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Abstract 

The fourth, missing example of an exactly solvable complex potential with VT sym- 
metry V{x) = [V(— x)]* defined on a bent contour and leading, at the real energies, to 
the Jacobi polynomial wave functions is found in a generalized Hulthen interaction. 
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Quantum mechanics is often forced to formulate its predictions numerically. Its ex- 
actly solvable models are scarce. Only their subclass in one dimension is broader 
and, in this sense, privileged and exceptional. It involves the harmonic oscillator and 
Morse potentials (with bound states expressible in terms of the Laguerre polynomi- 
als) as well as several models solvable in terms of the polynomials of Jacobi (cf., e.g., 
review |l| for more details). 

Within the framework of an alternative or modified quantum mechanics as pro- 
posed recently by Bessis |2| and by Bender et al |3], Q one works with the VT 
symmetric complex potentials V(x) = [V(— x)]*. The real bound state spectrum and 
even the solvability of several real potentials proves rather unexpectedly preserved 
after a complexification of this type. For example, within the Laguerre-related ex- 
actly solvable sub-family we may find both the complex harmonic oscillator || and 
a complexified Morse interaction ||. The same high degree of analogy between the 
real and complex forces is also observed for the regular solvable models based on the 
use of the Jacobi polynomials |J. 

A word of warning against unlimited optimism comes from singular interac- 
tions. In particular, for the Laguerre-related and phenomenologically most appealing 
Coulombic interaction the only available VT symmetrization proves merely partially 
solvable 0. In the present communication we intend to offer a partial remedy. We 
shall derive and describe a complete and exact solution for an appropriate VT sym- 
metric complexification of the singular phenomenological Hulthen potential which is 
known to mimick very well the shape of the Coulombic force in the vicinity of its 
singularity ||. 

In the first step let us recollect that in the spirit of the old Liouville's paper |J 
the change of the (real) coordinates (say, r < — ► £) in Schrodinger equation 

d 2 1 

--^ + w ( r ) X{r) = -K 2 x{r) (1) 
may sometimes mediate a transition between two different potentials. It is easy to 



show [10] that once we forget about boundary conditions one simply has to demand 
the existence of the invertible function r = ?"(£) and its few derivatives r'(£), r"(£), . . . 
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in order to get the explicit correspondence between the two bound state problems, 
viz., original eq. (P and the new Schrodinger equation 



d 2 



with the wave functions 



*(0 = x[r(0]/\MO 

generated by a new interaction at new energies, 
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(2) 



(3) 



V(0 -E= [r'(Of {W[r(0] + k 2 } + 
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(4) 



The mapping between the Morse and harmonic oscillators is one of the best known 
explicit illustrations of this rule since the necessary preservation of the correct phys- 



ical boundary conditions is very straightforward to check there [|TTJ . An appropriate 
VT symmetrized extension of this equivalence to the complex non-Hermitian cases 
is easy @. 

In the alternative, Jacobi-polynomial context the Liouvillean changes of variables 
have been applied systematically to all the Hermitian models (cf. Figure 5.1 in the 
review [I] or refs. |12| for a more detailed illustration). A similar thorough study is 
still missing for the VT symmetric models within the same subclass. 

In the present letter we shall try to fill the gap. For the sake of brevity we shall 
only restrict our attention to the VT symmetric initial eq. (|l|) with the Poschl- Teller 



potential studied and solved exactly in our recent preprint [13 
/5 2 - 1/4 a 2 -1/4 



W ( r ) = — T2 72— , r = x-ie, 

sinn r cosh r 

The normalizable VT symmetric solutions 



x G (— oo, oo) 



(5) 



X (r) = sinh T/3+1/2 r cosh ffQ+1/2 r P^'™) (cosh 2r) 

are proportional to the Jacobi polynomials at all the negative energies — k 2 < such 
that 

k = k^' t) = -era - t(3 - 2n - 1 > 0. 
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These bound states are numbered by n = 0,1,..., n^J and by the generalized 
parities a = ±1 and r = ±1. 

We may note that our initial VT symmetric model (p]) remains manifestly regular 
provided only that its constant downward shift of the coordinates r = rr x \ = x — is 
remains constrained to a finite interval, e G (0, 7r/2). In this sense our initial model 
(H) is closely similar to the shifted harmonic oscillator. At the same time, one still 
misses an analogue of a transition to its Morse-like final partner V(£) in eq (0). In a 
key step of its present construction let us first pick up the following specific change 
of the axis of coordinates, 

sinh r( x ) (£) = — ie 1 ^, £ = v — iu. (6) 

The main motivation of such a tentative assignment lies in the related shift and 
removal of the singularity (sitting at r = 0) to infinity (u — > +oo). In fact, one 
cannot proceed sufficiently easily in an opposite direction, i.e., from a choice of a 
realistic V(£) to a re-constructed r(£). This is due to the definition (f|) containing 
the third derivatives and, hence, too complicated to solve. 

We shall see below that we are quite lucky with our purely trial and error choice 
of eq. @. Firstly, we already clearly see that the real line of x becomes mapped 
upon a manifestly VT symmetric curve £ = v — iu in accordance with the compact 
and invertible trigonometric rules 

sinhxcose = e u sinv, cosh a; sine = e u cosv, 

i.e., in such a way that 

v = arctan (^f ) = v {x) G (^(-oo)^(oc)) = (-f+e,f-e), 
u = U(x) = \ In (sinh 2 x + sin 2 . 

This relationship is not too different from its Morse-harmonic predecessor studied in 
ref. ||. Our present path of £ is a very similar down-bent arch which starts in its 
left imaginary minus infinity, ends in its right imaginary minus infinity while its top 
lies at x — v — and — u = —uto) = lnl/ sine > 0. The top may move towards the 
singularity in a way mimicked by the diminishing shift e — > 0. Indeed, although the 
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singularity originally occurred at the finite value r — > 0, it has now been removed 
upwards, i.e., in the direction of — u — > +00. 

The first consequence of our particular change of variables (H) is that it does not 
change the asymptotics of the wave functions. As long as r'(£) = itanhr(£) the 
transition from eq. ([I]) to (|2|) introduces just an inessential phase factor in 
This implies that the normalizability (at a physical energy) as well as its violations 
(off the discrete spectrum) are both in a one-to-one correspondence. 

The explicit relation between the old and new energies and couplings is not too 
complicated. Patient computations reveal its closed form. With a bit of luck, the 
solution proves non-numerical. The new form of the potential and of its binding 
energies is derived by the mere insertion in eq. (0), 

At the imaginary £ and vanishing A = this interaction coincides with the Hulthen 
potential. 

In the new formula one has to notice the positivity of the energies. It is extremely 
interesting since the potential itself is asymptotically vanishing. One may immedi- 
ately recollect that a similar paradox has already been observed in a few other VT 
symmetric models || |14|, [15] where even an asymptotic decrease of the potential to 
minus infinity did not destroy a lower bound of the spectrum. 

The exact solvability of our modified Hulthen potential is not yet guaranteed at 
all. A critical point is that the new couplings depend on the old energies and, hence, 
on the discrete quantum numbers n, a and r in principle. This could induce an 
undesirable state-dependence into our new potential. Vice versa, the closed solvabil- 
ity of the constraint which forbids this state-dependence will be equivalent to the 
solvability at last. Fortunately, a smooth removal of the obstacle is possible by a 
transfer of the state-dependence (i.e., of the n— , a— and r— dependence) in 

A = A{a) = l-a 2 , C (= A + B) = k 2 - f3 2 

from C to (3. To this end, employing the known explicit form of k we may re-write 

C = C(a, r, n) = (a a + 2n + I) (a a + 2n + 1 + 2t(3). (8) 
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This formula is linear in r(3 and, hence, its easy inversion defines the desirable state- 
dependent quantity (3 = (3(a, r, n) as an elementary function of the constant C. The 
whole new energy spectrum acquires the closed form 



E — E(a, t, n) = A + B + - 



A + B 

aa + 2n + 1 — 



aa + In + 1 



(9) 



Our construction is complete. The range of the quantum numbers n, a and r remains 
the same as above. 

In the light of our new result we may now split the whole family of the exactly 
solvable VT symmetric models in the two distinct categories. The first one "lives" on 
the real line and may be represented or illustrated not only by the popular Laguerre- 
solvable harmonic oscillator |Hd] but also by our initial Poschl- Teller Jacobi-solvable 
force (H). 

The second category requires a narrow arch-shaped path of integration which all 
lies confined within a vertical strip. It contains at last both the Laguerre and Jacobi 
solutions. The former ones may be represented by the complex Morse model of ref. 
||. Our present new Hulthen example offers its Jacobi solvable counterpart. The 
scheme becomes, in a way, complete. 

The less formal difference between the two categories may be also sought in their 
immediate physical relevance. Applications of the former class may be facilitated by 
a limiting transition which is able to return them back on the usual real line. In 
contrast, the second category may rather find its most useful place in the methodical 



considerations concerning, e.g., field theories and parity breaking [I7|. Within the 
quantum mechanics itself an alternative approach to the second category might also 
parallel studies ITBI of the "smoothed" square wells in the non-Hermitian setting. 

In the conclusion let us recollect that the VT symmetry of a Hamiltonian replaces 
and, in a way, generalizes its usual hermiticity. This is the main reason why there 
existed a space for a new solvable model among the singular interactions. At the same 
time, an exactly solvable example with another, intermediate (i.e., hyperbola-shaped) 
arc of coordinates remains still to be discovered. Indeed, this type of a deformed 
contour has only been encountered in the quasi-solvable (i.e., partially numerical) 
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model of ref. [JH| ) and in the general unsolvable forces studied by several authors 
by means of the perturbative [19|, numerical [20] and WKB j|, [H[ approximation 
techniques. 
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